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ABSTRACT 

The paradigm that gamma-ray burst (GRB) firebaUs are the sources of the ultra-high energy cosmic 
rays (UHECRs) is being probed by neutrino observations. Very stringent bounds can be obtained from 
the cosmic ray (proton)-neutrino connection, assuming that the UHECRs escape as neutrons. In this 
study, we identify three different regimes as a function of the firebah parameters: the standard "one 
neutrino per cosmic ray" case, the opticaUy thick (to neutron escape) case, and the case where leakage 
of protons from the boundaries of the shells (direct escape) dominates. In the optically thick regime, 
photomeson production is very efficient, and more neutrinos will be emitted per cosmic ray than in 
the standard case, whereas in the direct escape-dominated regime, more cosmic rays than neutrinos 
will be emitted. We demonstrate that, for efficient proton acceleration, which is required to describe 
the observed UHECR spectrum, the standard case only applies to a very narrow region of the fireball 
parameter space. We illustrate with several observed examples that conclusions on the cosmic ray- 
neutrino connection will depend on the actual burst parameters. We show that the definition of the 
pion production efficiency currently used by the IceCube collaboration underestimates the neutrino 
production in the optically thick case. Finally, we point out that the direct escape component leads 
to a spectral break in the cosmic ray spectrum emitted from a single source. The resulting "two- 
component model" can be used to even more strongly pronounce the spectral features of the observed 
UHECR spectrum than the dip model. 

Subject headings: Gamma-ray burst: general. Methods: numerical. Neutrinos 



1. INTRODUCTION 

Cosmic ray observations tell us that particles with en- 
ergies higher than 10^^ eV hit the Earth, which are ex- 
pected to be of extragalactic origin. The search for the 
sources of these ultra- high energy cosmic rays (UHECRs) 
is therefore one of the main objectives in high-energy as- 
trophysics. It can be either performed directly, by cos- 
mic ray observations, or indirectly, by looking for the 
neutrinos accompanying the cosmic ray emission. So far, 
no evidence for a correlation between specific UHECR 
sources and cosmic ray measurements has been found. 

One class of potential UH ECR s o urces ar e gamma- 
ray bu rst (GRB) fireballs f see iPiranI ()2004[ ) or iMeszaroi 
()2006f ) for reviews), where the cosmic rays are ex- 
pected to be accelerated to the highes t energies by col- 
lision s with the interstellar medium (jRees &: Meszaro^ 
11992 ), or by internal collisions inside the ejected mate- 
rial ((Paczvnski fc Xu|[T99l IRees fc Meszarosl[l99l . The 



general fireball model describes a GRB as a catastrophic 
release of energy, during which matter of the order of 
a solar rest mass is ejected from a compact object. 
These ejecta are then first accelerated to ultra-relativistic 
speeds, and are then assumed to coast at constant ve- 
locity while expanding into interstellar space. It is as- 
sumed that the bursts can also lose energy via radia- 
tion at higher radii during the expansion of the fireball. 
The expansion ends when the ejecta hit the interstellar 
medium and are decelerated. The observations in sev- 



eral energy bands have shown that GRBs have several 
distinct phases of emission, with the prompt emission 
phase being the most energetic one. During this phase, 
the burst is mainly visible in gamma-rays, and the emis- 
sion is assumed to originate from the coasting phase of 
the fireball. The high variability and non-thermal prop- 
erties of the observed gamma-ray spectra give rise to the 
notion that the prompt emission might be due to the 
collision of internal shells, leading to Fermi shock ac- 
celeration of the charged particles such as electrons or 
protons (internal shock model). While recent observa- 
tions poi nt towards a heavier c omposition at the highest 
energies ^Abraham et al.ll20Tol ). we focus on protons as 
candidates for the UHECRs in this study, for which plau- 
sible models for the particle acceleration and emission 
from GRBs exist. Especially the idea of indirect escape 
of protons, via neutrons, from the shells has been a very 
popular extension of the internal shock model in the lit- 
erature. However, there are also several other alternative 
models for GRB emission which describe certain aspects 
of the observed emission well, such as mag netic reconnec- 
tion models ([ Lvutikov & Blandford 2003; iZhang fc YanI 
120111 : IZhang fc Kum ar 201 ^) or pho tosphcric emission 
models (IRees fc Me szaros 2 0051 iGiannios 20081: iMurase 
200 1 IWang fc Dail 2009: B eloborodovDIo'lOl: IGao eTal 



20TalLazzati et al.ll2013.) . 



Because of the high photon densities, it is expected 
that the protons accelerated in the colliding shocks dis- 
sipate energy into pion production. In the standard pic- 
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ture, this can be described by the A(1232)-resonance 



n + 7r+ i of ah cases 
p + I of aU cases 



(1) 



A substantial neutrino flux then originates from 7r+ de- 
cays via the decay chain 



+ i^e + i^fj. , 



(2) 



where 



i>T are produced in the ratio 1:2:0. On 



the other hand, the neutrons decay via 



(3) 



typically outside the source, which leads to a cosmic ray 
(proton) flux even if the protons themselves are mag- 
neti cally confined ("neutron model"), as discussed in, 
e.g., IMannheim et al.l (|2001[ ). Highly energetic gamma- 
rays originating from the tt^ decays are injected into 
the electromagnetic cascade, which leads to constraints 
from the Fermi-LAT diffus e GRB measurement; see, e.g., 
lAhlers et alJ (pOlOl) orO (|2012aD . 

The neutrino flux has b een predicted for th e stan - 
dard internal shock model in lWaxman fc Bahcalll ()1997[ ). 
assuming that GRBs are the sources of the UHECRs. 
If one assumes that the observed gamma-ray spec- 
trum represents the photon density within the source 
in the prompt phase, one can calculate the expected 
neutrino fluence from the observed gamma-ray flu- 
ence; corre spon ding analytical m ethods have b e en de - 
velo ped in iGue tta et al. (2004), iBecker et all (|2006D . 
and lAbbasi et al.l f|201Q) . Recently, the IceCube col- 
laboration has however str ongly constra i ned the neu - 
trino flux from GRBs, see lAbbasi et all (poTTl |2012| ). 
with the conclusion that these simple approaches are 
already severely constrained. Nonetheless, it is known 
that additional photomeson production processes some- 
what harden the neutrino spectra, and that the cooling 
of the secondaries and flavor mixin g change the spec- 
tral sh a pe and flavor compos ition, se elKashti fc WaxmanI 
(I2005D. iMurase fc NagataFj (2006), Lipari et al.l (I2007D . 
iHummer et al.l (|2010bl ). and ,Baerwald et al.l (|201lD . In 
addition, the predicted normalization is significantly re- 
duced if spectral effects on the pion production effi- 
ciency, the energy dependence of the mean free path 
of the protons, and the impact of the secondary cool- 
ing on the energy b u dge t are ta ken i nto account, se e 
IHiimmer et all (|20T1) . llJ (I2012b[) . and IHe et all (|20T1) . 
If one extrapolates a quasi-diffuse flux from a few GRBs, 
the low statistics of t he stacking sample wi ll lead to 
a systematical error ([Bacrwald_etaD l2012bD . Given 
the astro physical and systematical u ncertainties of the 
model in IWaxman fc Bahcallj (|1997D . the current neu- 
trino observations just start to enter the predicted neu- 
trino flux range, and the full-scale IceCube experiment 
should find neutrinos after ten years of operation if the 
baryonic loading of t he jets is as high as anticipated in 
IHummer et al.l (|2012[ ). Note, however, that if the prompt 
emission comes from larg er radii, no neutrinos may be 
found, see iHe et al.l (|2012h . Additionally, there have iDeen 
recent efforts to calculate the neutrino emissio n from 
GRBs in the dissipative photospheric models (Gao et al.l 



l2012f ), as well as efforts for model- independent calcula- 
tions (IZhang fc Kumadl20T2l) . 

As far as the direct connection between neutrinos and 
cosmic rays is concerned, equation ([l} suggests roughly 
one muon neutrino per cosmic ray after flavor mixing, 
which changes the Ve : : Vr from 1 : 2 : to 1 : 1 : 1 . 
Thi s hypothesis has be en tested in lAhlers et al.l ()2011[ ) 
and lAbbasi et al.l ()2012|) . with the conclusion that GRBs 
cannot be the sole source of the UHECR protons. In 
a more general framework, the authors of IKistler et al.l 
(2013:) conclude that the protons resulting from photo- 
pion processes are not sufficient to explain the cosmic-ray 
measurements. Therefore, we discuss the validity of the 
assumptions going into the "one neutrino per cosmic ray" 
paradigm, henceforth called "the standard case" : 

1. The protons are magnetically confined, and cosmic 
rays can only escape as neutrons. 

2. The protons interact only once, at most, and the 
produced neutrons can escape from the source 
(source optically thin to neutron escape). 

If one of these assumptions is violated, the consequences 
are obvious: protons "leaking" are not accompanied by 
neutrino production. On the other hand, multiple inter- 
actions will enhance the neutrino flux compared to the 
standard picture, while only neutrons from the bound- 
aries can escape. In this study, we will explore these two 
regimes in addition to the standard picture, and we will 
demonstrate that, for high proton acceleration efficien- 
cies, which are required to describe the observed UHECR 
spectrum, the standard case only occupies a very small 
region of the parameter space. 

This paper is organized as follows: We describe the im- 
plementation of the GRB fireball model for our purpose 
in Section |21 This section may be skipped if the reader 
is familiar with this or similar models. In Section [Sj we 
discuss the direct escape of cosmic rays from optically 
thin (to neutron escape) sources, and we comment on 
the additional effects of diffusion in Appendix [Xj On the 
other hand, we treat the optically thick (to neutron es- 
cape) regime in Sectional where we also comment on the 
pion production efficiency in that regime in Appendix [Bl 
In Section [5] we relate the important cases for the cosmic 
ray-neutrino connection to regions in parameter space, 
and we discuss specific (observed) examples in Section [HI 
Furthermore, we illustrate in Section |7| the impact of an 
additional escape component on the observed cosmic ray 
flux. Finally, we summarize in Section [8l 

2. IMPLEMENTATION OF THE GRB FIREBALL 
MODEL 

We use a simplified descri ption of the rela tivistically 
expanding fireball, based on iWaxmaiJ ()2003l) , to illus- 
trate our main points. Primed quantities refer to the 
shock rest frame (SRF), and unprimed quantities to the 
observer's or source (cosmologically co-moving) frame, 
which we clearly indicate. GRB observations exhibit a 
strong time variability over a scale ty (defined in the ob- 
server's frame), which can be related to a basic length 
scale ro = ctyl{\ -f z) in the source frame. We assume 
that the central engine of the GRB emits shells of thick- 
ness Ar ~ c-tv/{l -\- z) = tq in the source frame, since 
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causality implies that variations of the timescale can 
only propagate over a distance scale Ar. The time evolu- 
tion of the fireball can be divided into different zones. In 
the first zone, the shell gets accelerated, powered by the 
energy transfer from the thermal photons to the baryons 
in the shell. The Lorentz factor of the shell grows with 
the radius until a maximum value T is reached, which is, 
in principle, given by F = Etot/iMc^), where M is the 
total mass of baryons and -Etot, the total energy of the 
fireball. This transition is complete at a radius r « Fro- 
Here the second zone is considered to start: since the 
shell is accelerated to its maximal velocity, it coasts with 
constant F, while the expansion of the width of the shell 
itself is still negligible. However, when the shell reaches 
the radius r « F^ tq, the growth of the shell width can 
no longer be neglected, since 6Ar ~ r/F^ ~ tq. We will 
come back to this later when we discuss the effects of 
an expanding shell. For now, we assume that the shell 
width in the SRF is roughly given by 



Ar' ~ F c ■ 



1 + z 



at the indicated radius. At roughly the same radii, 
rc ^ 2T^ ro ^ 2T^ c 



1 + z' 



(4) 



(5) 



the collisions of the different shells start, based on the 
assumed fluctuations of the shells Lorentz factors of the 
order of AF/F ~ 1. External collisions with the interstel- 
lar medium can also lead to efficient proton acceleration, 
which we do not consider since the typical photon den- 
sities are orders of magnitude lower than in the internal 
collision zone. 

We focus on the description of the prompt phase, 
which is associated with the collisions of the shells. 
Since a relativistically expanding fireball may undergo 
different phases in its expansion, with varying param- 



eters, we des cribe t he phy s ics of one collision, f oUow- 
ing Baerwald gFaD (|2012bl) . iHiimmer etldl (I2012D. and 
IWinteri (120121) and cons i stent with iWaxman fc Bahcalll 



(|1997f) and lGuetta et"al1 (j2004[ ). If one assumes that the 
collisions occur at the same radius rc, as i t is implied in 
all of the st ate-of-the-art neutrino analyses (lAbb asi et al.l 
I20111I2012D . the total fiuences can be obtained by sum- 
ming over N ~ Tgo/t^ such collisions, where Tgg is the 
time during which 90% of the total energy is observed. 
Our shell-dependent approach has the advantage that the 
conventional results can be easily retrieved, and that in 
addition the relation to collision radius-dependent mod- 
els can be established. 

For the photohadronic interactions, the sec- 
ondary (such as pion) injection Q'{E') (in units of 
GeV^^ cm~^ s~^) can be computed from the proper 
photon N!^{e') and proton Np{E'p) densities (SRF, in 

units of GeV~^ cm~'^) as 



Q'{E') 



dEL 



EL 



^n;,{E'p) / cd£'7V;(e')i?(x,y). (6) 



Here, x = E' / E'^ is the fraction of energy going into 
the secondary particles, y = {E'pe')/mpC^ , and R{x,y) 
is a "response function" . If many interaction types are 



considered, this response function can be quite compli- 
cated. Nevertheless, if it is known from particle physics, 
equation ([6]) can be used to compute the secondary in- 
jection for arb i trary proton and photon spectra; see 
IHiimmer et al.l ()2010b( ). Note that the secondary injec- 
tion depends on the product normalization of the proton 
density Np{Ep) and the target photon density N'^[e'). 
Once the proper proton and photon densities (including 
the spectral shapes) are known, as well as the magnetic 
field -B', the secondary meson and neutron production is 
just a straightforward particle physics co nseq uence. We 
use the method from Hiimmcr et al. (20 10b( ). based on 
the physics of SOPHIA (Miicke et al.n200d[ rfor the com- 
putation of the photohadronic interactions. For the sec- 
ondary meson decays (includin g the helicity dependence 
of the muon decays), see, e.g.. iLipari et al.l (|2007 ) . The 
magnetic field effects an d fla vor mixing are included as in 
iBaerwald et"an (|2012bf ) and IHiimmer et al.l (|2012( ). Be- 
low, we will describe how to determine the relevant input 
N^iE'p), N^ie'), and B' from the observables. 

In equation ©, two types of spectra are present: the 
injection/ejection spectrum Q' and the steady spectrum 
N'. For a specific particle species, these are related to 
each other by a kinetic equation describing energy losses 
and escape. If the energy losses can be neglected, they 
are, for one species of particles, related by 



N'{E')^Q'iE')t',, 



(7) 



where t'^^^ is the escape time. For example, the observed 
gamma-ray spectrum can be obtained from Q'^ , whereas 
the spectrum relevant for the photohadronic interactions 
in equation ^ is N!^. Typically, one establishes a rela- 
tionship between observed gamma-ray fiuence and target 
photon density by implying that the gamma-rays escape 
over t^g^ = ijyjj ~ Ar' /c, which means that equation 
can be used. However, if the optical thickness to pair 
production or other processes is of order unity, this as- 
sumption does not apply, and the observed spectrum is 
not representative for the density in the source anymore. 

For 7V^(e'), a broken power law is normally assumed, 
parameterized as 



V ^7, break / 







pi < f' < f' 

''7, mill — ^ ^-7, break 

f' < f' < f' 

7, break — ^ '-7, max 

else 



(8) 



with e'^ — C'(keV) the break energy of the photon 
spectrum in the SRF. Typical values for the spectral in- 
dices are ss 1 and w 2. The minimal and the maxi- 
mal photon energies are chosen to be e!!y,niin = 0.2 eV and 
£7, max — 300 • |-,;.g3^j, in our calculations, if the highest 
energetic photons can escape (see below). These values 
are far enough away from the break energy t o have no vis- 
ible effect on the predicted n eutrino spectra ()Lipari et al.l 
120071: IBaerwald et al.|[20T2H ). though they can somewhat 
affect the neutron escape spectra; see Section [S] for a 
more detailed discussion. In addition, the energy parti- 
tion is hardly affected by the maximal photon energy for 
a spectral index > 2, since the energy in photons then 
depends only logarithmically on the maximal photon en- 
ergy, at most. Note that high-energy photons will not 
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be able to escape above the pair production threshold. 
In this case, we choose e^^max consistent with the pair 

production thresholdQ 

In case of the internal collisions, it is generally assumed 
that Fermi shock acceleration leads to a non-thermal par- 
ticle spectrum of the form 



exp 



E' 

p.max 



(9) 



with the spectral index ap sa 2. For the exponential cut- 
off, we choose k ~ 2 unless noted otherwise. The maxi- 
mal proton energy E'p .^^^ can be obtained by comparing 
the acceleration timescale to the dominant loss timescale 

^acc (-^p,max) min [^clyi;i , ^gyn (-^p.max) 7 (-^p,inax)] ' 

Here we assume that the acceleration time (in Gaussian 
cgs units) is given by 



E' 



TjceB' 



(11) 



with the elementary charge e ~ 4.803 • 10^^*^ Fr and rj 
the acceleration efficiency {rj is defined here so that large 
values mean efficient acceleration). It is generally as- 
sumed that the dominant loss mechanism is either the 
adiabatic loss due to the expansion of the shell or the 
synchrotron loss of the protons due to the magnetic fields 
present in the shells. We do not consider the adiabatic 
loss timescale explicitly, since we assume that it is of the 
same order as the dynamical timescale t'^yj^{E'). The 
synchrotron loss time is given by 



Cn(-^') - 



4ce4 E' 



(12) 



with the particle mass m being in erg, using the rela- 
tion 1 erg — 624.15 GeV. Moreover, the photohadronic 
timescale t'^^ is numerically computed from the interac- 
tion rate as given in iHlimmer et all (|2010b[ ). 

Let us now derive N!y{e') and Np{E'p) from the observ- 
ables. Frequently used observables are the (bolometri- 
cally corrected) gamma-ray fluence of a detected GRB, 
'5'boi (in units of ergcm"^), or the radiative flux (in 
units of GeVcm"^ s~^). Here we focus on a momentary 
picture of the fireball, described by (a possibly bolomet- 
rically corrected) i^-y, which leads to the isotropic equiv- 
alent energy per shell (or collision) 



«i, 47r d r 



(1 + ^) 



(13) 



^iso 



where is the luminosity distance. One has E[l 
Ef^^/T in the SRF, and i^,i,o = Ef^^ ■ (l + z)/t^. Assum- 
ing energy equipartition between photons and electrons. 



^ We use equation (6) from IWaxman fc BahcaHI 119971) to esti- 
mate that, applicable for the ^"■^-spectra (above the break), assum- 
ing that the photon spectrum extends to infinitely high energies. 
This is only a rough estimate, since gamma-rays may interact by 
additional processes. The impact on max however, typically 
small. 



the photons carry a fraction (fraction of energy in elec- 
trons) of the total energy E'fjj, j^^, and 



E\ 



ish 



— 1 psh 



(14) 



In order to compute the photon and proton densities in 
the SRF, it turns out to be useful to define an "isotropic 
volume" V^gQ = 47r Ar' cx F^ , where the latter relation- 
ship can be easily read off from Eqs. dH) and ([S]). Here 
T^gQ can be interpreted as the volume of the interaction 
region assuming isotropic emission by the source If the 
characteristics of all collisions are alike, S'boi — ^7 Tqq. 

Now one can determine the normalization of the pho- 
ton density in equation ([8]) and the proton density in 
equation ^ from 



e'7V'(£')de' = 



E: 



i/sh 



VL 



1 rp/sh 

e'^n;{e'^)ae'^= 



(15) 

Here /e is the ratio between energy in electrons and pro- 
tons (/~^ is the baryonic loading). Assuming that the 
magnetic field carries a fraction es of E[^, one has in 
addition 



or 



B' = \ 8tt 



ISO 

V 

ISO 



(16) 



After photohadronic interactions and weak decays, one 
obtains the injection spectrum of secondary neutrinos 
or neutrons Q' , which is to be translated into the ob- 
servable neutrino or neutron fluence J-^^ (in units of 
GeV~^ cm~^) per shell: 



F — f V' 



(1 + ^)^ 

47rdi 



E = 



1 



E' . (17) 



Note that for the sake of comparability, we show all neu- 
trino and cosmic ray fluences in the observer's frame, 
assuming that diffusion, pair production, and photo- 
hadronic interactions can be neglected. In Section [3 
however, we discuss the impact on UHECR observations 
including pair production and photohadronic losses. Ad- 
ditionally, we assume that neutrinos are subject to flavor 
mixing using the mixing angles 6x2 = 0.587, ^13 = 0.156, 
62?, = 0.670, and (5cp = 1-08 tt, taken from iFogli et al.l 
(201^ for the normal (mass) hierarchy. For reference, 
it will be also illustrative to show the cosmic ray proton 
fiuence if all protons were allowed to escape over t'^y^, 
which represents the maximal possible "leakage" from 
the source. This fiuence can be obtained from equa- 
tion m using Q; = N'Jt'^y^. 

There are three important features of our approach. 
First of all, we relate everything to the prompt phase, 
which is implied by using the flux during that phase in 
equation (|13p . The proper densities in equation ()15|) de- 
scribe the (steady) proton and photon densities in that 
phase. We do not specify the origin of the target pho- 
tons, such as synchrotron emission of co-accelerated elec- 
trons or inverse Compton scattering. Second, we con- 
sider emitted neutrino and cosmic ray fluences instead of 

^ Since both the energy and the volume of the source need to 
be, in principle, corrected by a beaming factor, this beaming factor 
cancels in the computation of the energy densities. 
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fluxes, which imphes that we do not need to resolve the 
time-dependence of the emissions. For instance, the cos- 
mic ray protons emitted with a different escape mecha- 
nism may not be emitted at the same time as the gamma- 
rays. And, third, we compute the fluences per shell, 
which may seem a bit peculiar, but has the advantage 
that our approach can describe dynamical changes of 
the fireball parameters over time, such as collisions at 
different radii. If all collisions are alike, as it is often as- 
sumed, one can easily obtain the result by summing over 
N ~ T^o/t^ such collisions. 

3. DIRECT ESCAPE FROM AN EXPANDING 
SHELL 

In this section, we discuss "direct escape" as the 
UHECR escape mechanism. This escape mechanism 
refers to the escape of protons without further inter- 
actions. Since protons (at lower energies) are magnet- 
ically confined, it is clear that only protons from the 
outer edges of the shells can escape directly. While this 
contribution may be generically expected to be small, 
we will demonstrate that it is an energy-dependent frac- 
tion of protons which can directly escape, and that the 
direct escape can dominate over the escape of neutrons 
produced in photohadronic interactions in regions of the 
parameter space. 

We describe the direct escape of protons from an ex- 
panding shell by the fraction of particles which can es- 
cape, relative to the densities in equation ()15|) . Assum- 
ing that the particles are isotropically distributed in the 
shell, the number of escaping particles is proportional to 
the escape volume. We assume that particles can escape 
from within a shell with thickness X[^^[p without inter- 
action, where AJ^fp refers to the mean free path. The 
fraction of escaping particles /esc present in the collision 
phase can be estimated from the ratio between escape 
volume and isotropic volume as 



v 



1 47r (r2 (r - Ar')2) a;, 



mfp 



A' 



mfp 



An 7-2 Ar' 



Ar' 



(18) 

Here it is taken into account that there are in fact two 
escape regions in each shell (inner and outer edge), and 
that only half of the particles along the mean free path 
propagate in the outwards direction. The proton flucnce 
from escape can then be computed from equation (|17p 
using 



N' 

dyn 



t': 



dyn 



Ar 



t' 



(19) 



cff.dii 



with an effective direct escape timescale ieffdir = 

^dyn^'^'/'^mfp- 

The relevant question is: what determines A^^fp of the 
protons? Since we consider GRBs, protons (or other 
charged particles) will interact with the magnetic field 
of the plasma. This is most often thought to trap the 
protons inside the shock due to magnetic confinement, 
while neutral particles, such as neutrons, are considered 
to escape. When the interactions with the magnetic field 
are treated as scattering, the Larmor radius 



R'^^~ 33.3cm- ' ^ ' 
e B' 



lOHJ 
B' 



(20) 



can be used to estimate the mean free path. That means 
that effectively everything within i?^ of the edges will 
eventually escape from the shells, which we call "direct 
escape" . Using the analogy to scattering, it corresponds 
to the region where the particles can escape without be- 
ing scattered at all, and without having lost energy. Now, 
R'j^ is proportional to and i?^ = ct'^^^ for 77 = 1 (see 
equation ([Tl]) V Therefore, it is easy to see that all pro- 
tons will directly escape at the highest energy if 77 = 1 
and the maximal proton energy is limited by the dynami- 
cal timescale, where ct'^^^ = Ar'. This is the region where 
direct escape of UHECR dominates. For lower acceler- 
ation efficiencies, or synchrotron- or photohadronic-loss 
limited maximal proton energies, the direct escape com- 
ponent will be smaller, and the question of what kind of 
escape dominates will be more complicated. 

As long as the proton and photon densities are not too 
high, the mean free path of protons is given by 



(21) 



This can be used in equation (jlSp . where the first term 
ensures that /esc < 1- Note however that this only ap- 
plies as long as the source is optically thin. We will 
discuss the optically thick regime, in which the photo- 
hadronic interactions dominate, in the next section in 
greater detail, since it turns out that it is different from 
the direct-escape dominated region. 

After having identified the different processes which 
limit the mean free path of the protons, we still need to 
consider how the expansion of the shell affects the scat- 
tering treatment of the particles. Beyond the collision 
zone, the shell width expands as Ar' cx r((x t). As a 
consequence, the volume of the shell grows as V' fx . 
This expansion is generally taken to be adiabatic since at 
this stage the photons are assumed to no longer couple to 
the electrons (or protons) . The general relation between 
energy and volume of a gas is given by 



E' oc y'-(^-i) oc r-3(7-i) 



(22) 



for the adiabatic index 7. Since we know the scaling of 
the volume as well as the total energy, we can also derive 
the scaling of the magnetic field inside the plasma as 




cx r 2 



(23) 



If we assume that the energy of a single particle in the 
plasma scales in the same way as the total energy of 
the plasma (see Eq. ([22l) ). the Larmor radius i?^ of the 
particles changes with time (or radius r) as 



(24) 



The usual assumption for a relativistic gas (or plasma) 
is 7 ~ 4/3, which leads to i?^ cx r. This scaling of the 
Larmor radius has an interesting consequence: the ra- 
tio between i?^ and the width of the shell Ar' cx r is 
constant. By computing neutrino and proton fluences 
(instead of fiuxes), we do not need to identify when the 
particles actually escape, we just compute the fraction of 
escaping particles /esc- This means that we can evaluate 
equation (jl8p in the collision phase using r ~ rc with 
the corresponding proton density. Note that other pos- 
sible values for 7 are 7 = 5/3 (mono- atomic classic gas). 
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7/5 (diatomic classic gas), or 9/7 (three-atomic gas with 
non-static bindings). Prom equation p4|. one can easily 
see that 7 = 4/3 represents a special case that simplifies 
the scaling and just gives stability between losses and 
escape. For 7 > 4/3, the adiabatic energy loss is faster 
than the escape, and the particles are more strongly con- 
fined for larger radii. In that case, our computations may 
over-estimate the direct escape, and most particles will 
directly escape early. For 7 < 4/3, the Larmor radius 
increases faster than the particles lose energy, and the 
particles are getting less confined at larger radii. With 
larger radii, particles with lower and lower energies can 
finally escape. In that case, especially lower energy pro- 
tons may escape in later phases. 

Two examples of direct escape-dominated bursts can 
be found in the different columns of Figure [T] for the pa- 
rameter sets given in the plot. The upper row shows 
the particle fluences, where "initial p" stands for the 
case if all protons were able to escape over the dy- 
namical timescale, "CR from n" represents the cosmic 
rays through neutron escape from photohadronic inter- 
actions, "direct escaping p" for direct proton escape, and 
+ ^11 for the muon neutrino fluence including flavor 
mixingO In the lower panels, the acceleration timescale 
and the considered energy loss/escape timescales are 
shown, where the direct escape effective timescale is de- 
fined in equation (|19p . It can be read off from these 
panels that in both cases the dynamical timescale lim- 
its the maximal proton energy. In the upper left panel, 
the direct escape clearly dominates. Since the accelera- 
tion efficiency ?? = 1, practically all protons escape at the 
highest energy, which is where the dynamical timescale, 
acceleration timescale, and direct escape timescale meet. 
However, in that case hardly any neutrinos are pro- 
duced due to the low photohadronic interaction rate. 
The upper right panel represents the typical case for the 
optically thin (to neutron escape) source, where both 
substantial neutrino and neutron fluxes are produced. 
The additional component from direct escape still dom- 
inates at the highest energies, while for energies below 
10^'^ GeV, the neutron flux dominates the cosmic ray 
production. This can be also read off from the corre- 
sponding timescales in the lower panel. Note that direct 
escape strongly depends on the acceleration efficiency: if 
77 ^ 1, the Larmor radius will be much smaller than Ar' 
at the maximal proton energy, and the direct component 
becomes suppressed. 

Note that the results in this section can be only inter- 
preted as rough estimates, and there may be additional 
escape components compared to the ones discussed here, 
e.g., diffusion may play a role. We discuss the possible 
impact of diffusion in Appendix [Xj where we demon- 
strate that it does not affect our qualitative conclusions. 
However, we also point out that a dedicated treatment 
of diffusion requires a mo del- dependent solution of the 
transport equations, which goes beyond the scope of this 
study, whereas direct escape can be regarded as a guar- 
anteed contribution to the cosmic ray injection. 



^ Note that for the sake of comparability, the CR spectra shown 
here are "at the observer" assuming that the CRs receive the same 
boost and losses as the neutrinos. In particular, the spectra here 
are without any losses during propagation apart from the adiabatic 
losses due to the cosmic expansion. 



4. OPTICALLY THICK (TO PROTON AND 
NEUTRON ESCAPE) CASE 

We define the optical thickness to neutron escape as 



_ >7 



f 



dyn 



(25) 



at the maximal proton energy. Thus, if t,i > 1, neu- 
trons at the maximal proton energy will rather inter- 
act than escape, and will therefore be confined. Since 
increases with energy (see, e.g., Figure[Tl lower pan- 
els), this optical thickness is typically at its maximum 
at the maximal proton energy. That is, it applies to the 
UHECR part of the emission, whereas neutrons at lower 
energies may escape more easily. Of course, not only the 
neutrons will interact rather than leave the region, but 
also the protons. 

For a numerical description, we follow the same ap- 
proach as in the previous section. We assume that only 
a fraction of the produced neutrons can directly escape, 
i.e., we apply the same mechanism to (direct) neutron 
escape and multiply the neutron injection Q'^ by equa- 
tion (1181) . The mean free paths for protons and neutrons, 
respectively, are in the optically thick regime given by 

A;,„,fp(i?')-nnn [/\r' ,R'^{E'),ct'^^{E')\ , 

A;,n.fp(i^') = niin[Ar',ct;^(£;')] . (26) 

We show two examples for the optically thick case in 
the columns of Figure[21 In both cases (see lower panels), 
the photohadronic interaction rate exceeds the dynami- 
cal escape rate by a factor of t„ at the highest energies. 
It also limits the maximal proton energy in both cases. 
The neutron production is therefore very efficient; see up- 
per panels. However, only the neutrons from the edges 
can escape, which implies that the dashed curves (corre- 
sponding to the escape of all protons over the dynamical 
timescale) cannot be exceeded, and a level of about 50% 
of the dashed curves is reached (since the baryon system 
contains about 50% protons and neutrons each in the 
optically thick limit). In contrast, the neutrinos from in- 
teractions everywhere within the shell can escape, which 
means that the neutrino fluence becomes relatively en- 
hanced, and the "one neutrino per cosmic ray" paradigm 
does not hold anymore. This can be especially seen in 
the upper right panel. Note however that the neutrino 
fluence is typically lower than the cosmic ray fluence, be- 
cause the neutrons obtain a higher fraction of energy in 
the interaction. In none of the discussed cases the di- 
rect escape of protons substantially contributes, which is 
characteristic for the optically thick regime. 

Our approach has several limitations. First of all, one 
may argue that the neutrons trapped by photohadronic 
interactions may escape later in a relativistically expand- 
ing fireball. Indeed, since the photon density drops as 
r^^ , the injection of neutrons ceases and the confine- 
ment by photohadronic interactions will come to an end 
at a certain radius (the "neutronsphere" ) , and all re- 
maining neutrons may escape. The level of the neu- 
tron density is at the level of the proton density at 
the end of the confinement phase, which means that 
a neutron fluence up to the level of the dashed curves 
in Figure [2] might be reached. Therefore, we expect 
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Fig. 1. — Particle fluences per shell (upper row) and inverse timescales of different processes (in SRF, lower row) as a function of E in 
the observer's frame. The figure shows two different parameter values for -L-|,^iso in the different columns, where the other burst parameters 



are fixed to F = 300, tv = 0.01 s, 77 = 1, ee/fB = li /e = 0.1, a~f = 1, = 2, e'^ = 1 keV, and z = 2. Both examples are for the optically 
thin (to neutron escape) case, where the optical thickness t„ is given in the panels. For the cosmic rays, only adiabatic energy losses are 
taken into account for the propagation. See main text for details. 



that our "CR from n" curves in Figure [5] represent a 
more sophisticated time-dependent calculation within a 
factor of two. Furthermore, the assumed energy parti- 
tion fractions may be different than the ones assumed 
in the optically thick regime. W e also do not consider 
the effects of muon r e-acceleration (iKoers fc Wiiersll2007l: 
iMurase et al.ll2012t iKlein et al.ll2012l). or i nteractions of 
pious and kaons ( Kachelriess et al.l I2008D . which have 
however much smaller interaction rates than the protons. 

There are several subtleties in the optically thick case, 
which are best illustrated with the pion production effi- 
ciency relevant for neutrino production. As we show 
in detail analytically in Appendix [B] these lead to an un- 
derestimation of the neutrino production in the optically 
thick case if the curren t IceCub e method for the compu- 
tation of /tt from Ab basi et al.l (2010) is used, which is 
the found a,tion for all state-of-t he-art GRB stacking anal- 
yses, as in I Abbasi et al.l ()2012D . In fact, it turns out that 
the original formula for from lGuetta et all ()2004|) also 



applies to the optically thick case if the energy partition 
is defined with respect to the particle densities within 
the source, even though it was not derived for that limit. 
Current state-of-the- art numerical predictions, such as 
iHiimmer et al.l ()2012| ). take this into account automati- 
cally. 

5. PARAMETER SPACE STUDY OF THE COSMIC 
RAY-NEUTRINO CONNECTION 

For the cosmic ray-neutrino connection, we identify 
three different regimes: 

Optically thin to neutron escape regime: This is 
the usual scenario discussed in the literature: the 
cosmic rays are produced as neutrons and can es- 
cape the source ("neutron model"). Additional 
escape components are negligible, and the "one 
(muon) neutrino per cosmic ray" paradigm applies. 

Direct escape regime: Here the cosmic rays from di- 
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Fig. 2. — Particle fluences per shell (upper row) and inverse timescales of different processes (in SRF, lower row) as a function of E in 
the observer's frame. The figure shows two different parameter values for L-y^isa in the different columns, where the other burst parameters 
are fixed to the same values as in Figure [l] Both examples are for the optically thick (to neutron escape) case, where the optical thickness 
r„ is given in the panels. See main text for details. 



rect escape dominate at least at the highest energy. 
Since the neutron production by photohadronic 
processes is sub-dominant, the one neutrino per 
cosmic ray relationship does not hold, and more 
cosmic rays than neutrinos will be produced. See 
Section [31 

Optically thick to neutron escape regime: Here 
the protons and neutrons interact multiple times, 
and only protons and neutrons on the outer edges 
of the shells can (directly) escape. The neutrinos, 
however, can escape from everywhere within the 
shell, which leads to more neutrinos per cosmic 
ray than in the optically thin case. See Section S) 

In this section, we discuss the type of cosmic ray source as 
a function of the GRB parameters. We distinguish the 
dominant effect by using figures such as Figure [T] and 
Figured) if the fluence maximum in the spectrum comes 
from directly escaping protons, we assign the direct es- 
cape category, otherwise the optically thin regime. The 



optically thin and thick cases are distinguished by the 
optical thickness t„, as defined in equation (^51 . being 
smaller or larger than one, respectively. It turns out that 
either of these three categories can be uniquely assigned 
(neglecting minor overlap). 

We show a (numerical) parameter space scan in Fig- 
ure [3l where we always plot L-y,iso on the horizontal axis. 
Let us focus on moderate acceleration efficiencies 77 = 0.1 
first (lower row), which clearly exhibit all three regimes. 
The optically thin case can be found close to the often 
used standard parameter values. In this case, the di- 
rect escape contribution cannot reach the same level as 
shown in Figure [T] since the maximal achievable proton 
energy is lower, and therefore the Larmor radius at the 
maximal energy cannot reach the shell width. Therefore, 
the escape of neutrons produced in photohadronic inter- 
actions dominates. On the other hand, in the light red 
(gray) regions, r„ > 1, which means that the neutrino 
production is enhanced. In the blue (darker) regions, the 
direct escape component exceeds the neutron escape at 
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Fig. 3. — Parameter space scan of the GRB parameters for classification of regimes. The upper row corresponds to the acceleration 
efficiency r] = 1, the lower row to r} = 0.1. For the standard parameter values, see caption of Figure [T] Dashed curves mark the limit 
between optically thin and thick regime s. Th e thick dashed curve represents the analytical result if the photohadronic interaction rate based 
on equation (3) of Waxnian & Bahcal| II1997I ) is used. In the dark-shaded regions "LAT invisible" , gamma-rays above 30 MeV cannot leave 
the source anymore due to pair production. That is, sources left of these regions are in principle visible in the full Fermi-GBM range and 
may be even observable in LAT, whereas sources within these regions will not exhibit emission into the LAT range. The dots correspond 
to the parameter sets chosen in Figure [T] and Figure [2] 



the highest energies. For efficient acceleration, 77 = 1, see 
upper row. Here the opticahy thin to neutron escape re- 
gion almost vanishes, which is due to an enhancement of 
the direct escape. Remember that for 77 = 1, all protons 
can directly escape at the highest energy if the maximal 
energy is limited by the dynamical timescale. To summa- 
rize, for efficient proton acceleration, the standard case 
(one neutrino per cosmic ray) only applies in a very small 
region of the parameter space, and either fewer or more 
neutrinos per cosmic ray are produced, depending on the 
parameters. 

In order to better understand the relationship to ear- 
lier works, consi der the A-resonance parameterization 
proposed bv iWaxman &: Bahcalll (|1997f ) (see Eq. (3), in- 
creased by a factor of two because we consider the photo- 
hadronic interaction rate, not the pion production rate) . 
We show in Figure[3]the separator between optically thin 
and thick regimes as thick dashed curves for this cross 
section. Obviously, in all cases the optically thin region 
reappears and gets enlarged. One reason is that the full 
numerical computation includes high-energy processes, 
such as multi-pion production, which enhance the inter- 
action rate at high energies by a factor of a few, and 
which is not included in the shown A-resonance approx- 
imation. As a result, photohadronic processes become 
more important in the numerical result. 



Let us now relate the parameter space to Fermi-luAT 
observations. We hence show in Figure [3] the "LAT in- 
visible" regions, where gamma-rays above 30 MeV can- 
not leave the source anymore because they exceed the 
pair production threshold. That is, sources left of these 
regions should be visible in the full Fermi-GBM range 
and may be even observable in LAT, while sources in the 
"LAT invisible" regions will not exhibit emission into 
the LAT range. One can clearly see that the parame- 
ter set corresponding to the rightmost dot, associated to 
the right column of Figure [21 cannot be seen in LAT. 
In addition, even though the optical thicknesses t„ (for 
neutrons) and r^-y (for photons) are roughly proportional, 
they still have slightly different parameter dependencies!! 
Note that, in practice, the LA T emission lasts longer 
(see, e.g., iLange fc Pohll (|2013[ )'). which means that it 
may come from larger emission radii. That is, in a re- 
alistic time-dependent model, the fireball may follow a 
trajectory in the considered parameter space, and may 

This is mainly due to one important difference in the calcu- 
lations, namely which frame is considered to be relevant. The 
maximal proton/neutron energy is calculated in the SRF during 
our calculation, however the calculation of Tjj is done for an ob- 
served photon energy. Hence there is a difference in redshift z and 
Lorentz factor F. Moreover, the break energy is important for the 
calculation of the (energy) densities, but it is not relevant for the 
optical thickness of the observed gamma-rays, as this calculation 
refers to photons far above the break. 
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actually visit more than one regime during the burst du- 
ration. It remains to be seen what the consequences of 
such a time-dependent model on the neutrino fluxes will 
be. 

Since Figure |3] only shows sections through the param- 
eter space, it is instructive to at least have approximate 
analytical expressions for the different regimes. The in- 
terface between the direct escape and optically thin re- 
gions can be obtained from the maximal proton energy: if 
it is dominated by the dynamical timescale, direct escape 
will dominate at the highest energies since all protons 
can escape (for rj = 1); if it is dominated by synchrotron 
losses, neutron escape will take over. Taking into account 
?7, one can estimate that 



i^J:L°o*< 3.6- 10^1 ergs- 



ty 



10 



2.5 



0.01 s 



l + z 



(27) 



limits the direct escape dominated region. The interface 
between the optically thin and thick regimes (analytical 
thick dashed curves) can be obtained from t„ ~ 1 in 
equation (p5)) . using the analytical expression for t^"^ 
from iWaxman fc Bahcalll ()1997D . The optically thick 
regime is then estimated as 



L;%f^-^>l.l-10''ergs-'- ^^^^^ 



0.01 s 



7, break 



keV / 



(28) 



These formulas allow to estimate how a specific burst can 
be classified for arbitrary parameters. There are, how- 
ever, some limitations. First of all, equation (j28p under- 
estimates the photohadronic interactions, as discussed 
above. And, second, some of the (numerical) parame- 
ter dependencies in equation ([27]) cannot be reproduced 
within these assumptions. 

6. APPLICATION TO SPECIFIC GRBS 

In the previous sections, we have discussed the depen- 
dence of the cosmic ray escape on the theoretical pa- 
rameters, such as Lj/iso- However, actually the gamma- 
ray flux or fluence are the observables, and Lj^iso is just 
a function of these observables. In addition, for many 
bursts, the other necessary ingredients, such as redshift 
and time variability, have been measured. We there- 
fore study in this section the different cosmic ray escape 
mechanisms for specific bursts; see Table [TJ In this case, 
we use the bolometric fiuence as observable, and assume 
that it is obtained from N = Tgo/ty identi cal collisions, 
in consisten cy wit h the approach e s used in lAbbasi et al.l 
((20Tll[20Tl . and lHiimmer"eraII ([20T2l) . 

The standard burst "SB" has been inspir e d to produc e 
a spectrum similar to IWaxman fc Bahcalll (|1997l Il999f ). 
The other three bursts have been actually observed, and 
their properties can be taken from the literature (see ta- 
ble caption). GRB080916C is one of the brightest bursts 
ever seen, although at a large redshift, and one of the 
best studied Fermi-hAT bursts. The gamma-ray spec- 
trum of GRB090902B has a very steep photon spectral 



cutoff and a smaller redshift, although P is very high as 
well. GRB091024 can be regarded as a typi cal example 
repre sentative of many Fermi-GBM bursts ()Nava et al.l 
1201 Ih . except for the long duration. Note that the neu- 
trino spectra from these GRBs have been also discussed 
in Winter (2012) and Bacrwald et aL, C2012a). 

We show in Figure 2] the expected cosmic ray and 
neutrino fluences for the GRBs listed in Table [T] (thick 
curves). In each panel, we also give the optical thick- 
ness to neutron escape and photon escape (for 30MeV 
gamma-rays). Note that we use an acceleration efficiency 
7] = 0.1 for these simulations; a higher acceleration ef- 
ficiency would pronounce the direct escape component. 
The burst SB is a typical example for an optically thin 
burst with a substantial amount of neutrino production. 
The contribution of the direct escape component depends 
on the proton acceleration efficiency, and in this case it 
is suppressed. Therefore, the relationship one neutrino 
per cosmic ray holds. On the other hand, the two high-P 
Ferrm-hAT bursts GRB080916C and GRB090902B ex- 
hibit a clear direct escape domination, which is a feature 
of the large P; see Figure [3] (second row, second col- 
umn). It is clear that GRBs will be only observed in 
LAT if they are left of the dark-shaded regions in Fig- 
ure [31 An even larger value of the photon energy, such 
as 100 MeV, will extend these regions further to the left, 
which means that the LAT-associated parameter space 
region tends to overlap the direct escape regime. There- 
fore, LAT-observed GRBs tend to directly emit UHECR 
protons at the highest energies, at least in the LAT emis- 
sion phase. The situation is different for GRB091024 in 
the lower right panel of Figure SI which is on the edge of 
the optically thick regime. For this burst, again cosmic 
rays from escaped neutrons dominate, and the neutrino 
production follows the standard assumption. 

In order to illustrate the impact of the minimal and 
maximal photon energies in equation ([5]), we illustrate 
the impact of a "bolometric correction" as shaded ar- 
eas in Figure |4l This bolometric correction takes into 
account that the gamma-ray fiuence has only been ob- 
served in a certain energy range, whereas it cannot be ex- 
cluded that lower and higher energy photons are present 
in the source as well, either because they are outside the 
detection energyrange, or because they cannot escape 
from the sourceQ The proton density is then calculated 
from energy partition using the extended energy range 
according to equation (|15|) . The bolometric correction 
increases it in all cases due to photons not accounted for 
in the observation, and therefore the normalization of 
the spectra increases, including that of the "initial pro- 
tons" . In neither of the cases, the neutrino spectrum is 
very much affected by this bolometric correction, apart 
from the normalization change0 The extension of the 

^ For this correction, we fix the gamma-ray spectrum in the 
observed energy range from the observed fluence, and then hn- 
early extrapolate the spectrum (on a double log plot) to the range 
between 0.2 eV (SRF) and 100 MeV (observer's frame). This ex- 
tended range is motivated by the fact that high energy protons 
then always find sufficiently many low energy photons as inter- 
action partners, and low energy protons find enough high energy 
photons. The 100 MeV are chosen in the observer's frame since 
they correspond to a typical Fermi-hAT energy. 

^ This correction depends mostly on the upper spectral index 
of the photon spectrum. It is e'~^ in the upper left panel, which 
leads to a logarithmic dependence on the maximal proton energy, 
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TABLE 1 

Properties of four bur s ts dis c ussed in Sect i on [6l see Baerw.-vl d et al.I (I2011) for SB ("Standa rd Burst", similar to 
IWaxman fc BahcallI (IT997I . [T999D '). INava kFal] 1I2011I) and.Greiner et al.I (l2009i') for G RB080916C. INava et al.I 1I20T1I) and 
IAbdo et AL.I H2009I) FOR GRB090902B. and INava et al.I H2011I) and IGruber et alJ (|2011I) for GRB091024. The luminosity is 
calculated with L-yjso = 4:TTdj^ -SBOL/Tgo, with 5bol the fluence in the (bolometrically adjusted) energy range 

IkeV- IOMeV. Adopted from IHummerI II2013I) . 



photon energy range, however, hardly affects the neu- 
trino spectral shape. On the other hand, the extension 
of the photon spectrum to higher energies has a signif- 
icant effect on the neutron spectra at low energies. In 
the lower right panel, the bolometric correction even 
leads to a lower maximal proton energy, which is be- 
cause photohadronic energy losses take over to limit the 
maximal proton energy. Additionally, we checked that 
the gamma-rays from Tr'^-decays produced by the three 
actually measured bursts are below the LAT bounds or 
observations for these bursts. 

In summary, Fermi-LAT observed GRBs seem to a 
have a strong direct escape component of cosmic rays, 
at least in that emission phase, whereas the direct es- 
cape contribution of typical GBM bursts depends on the 
proton acceleration efficiency. Therefore, conclusions on 
the cosmic ray-neutrino connection will depend on the 
actual burst sample including the specific parameters of 
the fireballs, and the time evolution of the fireball proper- 
ties, which we did not discuss in this section. A possible 
bolometric correction beyond the observed energy ranges 
of the gamma-rays typically has a small effect, as long as 

7. IMPACT ON UHECR OBSERVATIONS 

So far, we have discussed the cosmic ray and neutrino 
emission from a single source in terms of its fluence. How- 
ever, especially for cosmic rays, single sources may not be 
resolvable, and the injection as a function of redshift de- 
termines the observed spectrum and cosmogenic neutrino 
flux. There are four key ingredients which determine 
the observed cosmic ray spectrum at the highest ener- 
gies: the injection spectral index in the co-moving frame 
(corresponding to our ap), the maximal particle energy 
(our -Ep.max), the evolution of the density of sources with 
redshift, and the chemical composition of the UHECRs. 
Here wc focus on protons as cosmic rays and a,ssurn e the 
standard GRB evolution ( Hopkins fc Beacoml ()2006D star 
formation rate including iKistler et al.l correction 
for GRBs), which leaves ap and Ep^^g.^ as free param- 
eters. For a recent review on c osmic ray transport and 
models, seelAloi sio et al.l (j20 1^ and references therein. 

In order to describe the UHECR spectrum, it is a nec- 
essary condition that high enough proton energies can be 
obtained in the co-moving (source) frame. We show in 

and steeper in the other cases, which leads to a (stronger) power 
law suppression. 



Figure [5] the maximal proton energy log2o(-E'p.inax/GeV) 
(contours) in the co-moving frame as a function of the 
GRB parameters, similar to Figured The upper row is 
shown for an acceleration efficiency rj — 1\ the lower row, 
for rj = 0.1. One can clearly see the symmetry around 
the transition curve between the optically thin (or direct 
escape dominated) and optically thick cases, where the 
proton energy is limited by photohadronic losses in the 
latter case. The highest proton energies can be obtained 
along this transition curves, but the overall dependence 
of the maximal proton energy on the model parameters 
is relatively weak. The main impact comes from the ac- 
celeration efficiency (compare upper with lower row) : for 
77 = 1, about an order of magnitude higher proton ener- 
gies can be achieved than for rj = 0.1. 

For E > 10^'^ GeV, the mean free path of the protons is 
of the order of 1 Gpc {z ~ 0.25), while for E > 10" GeV, 
it is of the order of 100 Mpc (z ~ 0.024); the correspond- 
ing parameter spaces are illustrated by arrows in the 
L^y^iso-z panels of Figure [5j Therefore, for the UHECR 
spectrum, only low redshifts are relevant, while for the 
cosmogenic neutri no flux, the ful l redshift range can con- 
tribute (see, e.g., iKotera et alJ (|2010l )). and especially 
the very highly energetic protons will lead to a substan- 
tial neutrino production due to energy losses on the cos- 
mic microwave and infrared backgrounds. From the red- 
shift dependencies in Figure [5] it is clear that for 77 = 1 
(upper row, left plot), all bursts for the chosen parame- 
ters can produce the UHECRs above about 10^° GeV; 
and only bursts between L^^so — 6 • 10^" ergs~^ and 
]^q53 above about 10^^ GeV (illustrated by ar- 

rows). For ?7 = 0.1 (lower row, left plot), the range 
6 • lO^'^ ergs^^ < L^^so ^ lO^'^ergs"^ can produce the 
UHECRs above about 10^° GeV, and energies as high 
as 10^^ GeV are difficult to reach. Note, however, that 
these parameter ranges are found for the other parame- 
ters fixed to their standard values. Nevertheless, if GRBs 
are the sources of the UHECRs, this discussion illustrates 
that either the acceleration must be very efficient, or only 
GRBs from a very narrow parameter space region con- 
tribute at the very highest energies. 

In the following, we adopt a very pragmatic point of 
view to show the impact of the direct escape on the ob- 
served cosmic ray spectrum: we choose a reasonable pa- 
rameter set for a GRB with a high enough proton en- 
ergy from a high acceleration efficiency rj ~ 1 which 
exhibits the direct escape component at the highest en- 
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Fig. 4. — The expected cosmic ray and neutrino fluences for the GRBs listed in Table [T] The maximal proton energy and the optical 
thickness are given/marked in the different panels. The thick curves correspond to the photon fiuence and energy range given in the caption 
of Table[T] whereas the shadings represent a bolometric correction beyond this range; see main text. In each panel, we also give the optical 
thickness to neutron and photon escape (the latter for 30MeV gamma-rays). Note that we use an acceleration efficiency rj = 0.1 for these 
simulations. 



ergies, as well as the standard neutron escape compo- 
nent at lower energies, such as in Figure (TJ upper right 
panel. Then we assume that all bursts are alike in the 
CO- moving frame, i.e., we use a universal cosmic ray in- 
jection function based on this burst and following the 
standard GRB evolution, with an arbitrary normaliza- 
tion. The burst parameters are given in the caption of 
Figure [51 and the maximal proton energy is found to be 
1.9- 10^"'^ GeV, with which the protons are injected into 
the interstellar medium. For the cosmic ray transport, 
we use a deterministic kinetic equation sol ver assuming 
continuous e nergy losses, which is based on lAhlers et al.l 
()2009l I2010D . After the cosmic ray propagation, the nor- 
malization of the cosmic spectrum is chosen to reproduce 
the HiRes data (A bbasi et al...2008,) . Note that the spe- 
cific weight of the direct versus neutron escape depends 
on the chosen parameter set, and that a more refined 
treatment of the leakage may lead to adjustments. How- 
ever, the idea here is that the presence of an additional 



escape component can lead to a spectral break in the 
cosmic ray spectrum from one source, as it is visible in 
Figure [TJ upper right panel. In addition, note that in 
principle information on the cosmic energy budget can 
be obtained from this procedure. Since a detailed dis- 
cussion goes beyond the scope of this study, it will be 
performed elsewhere. 

Now the only remaining free parameter is ap. In that 
context, two often used approaches for the reproduction 
of the UHECR shape are: 



The "dip model": using an injection index at the 
source ap > 2.5 (depending on the source evolu- 
tion). It can reproduce the dip and the ankle very 
well, i.e., the region above about 10^ GeV, but the 
large ap is difficult to motivate from Fermi shock 
acceleration. Note that we for now neglect the 
problems with the e nergy budg et in this model, see 
iMurase et all (|2008[ ) or .Berezinskv et al.l (|2006|) . 



13 




10™ lO" lO" lO" 10™ lO" lO" lO" 10™ lO" lO'' 10- 

L^iso/erg-i"' L^i^Jerg-s'^ L^i^^/erg-j"' 




Fig. 5. — Maximal proton energy logjQ(Ep,max/GeV) (contours) in the source (cosmologically co-moving) frame as a function of tlie GRB 
parameters. The upper row is for rj = I; the lower row, for rj = 0.1. Same regions and parameters as in Figure|3] 

10^ 



'l 10^ 

o 



10' 



8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 




Fig. 6. — Example of observed cosmic ray (proton) spectrum for different assumptions at the source and different spectral injection 
indices ctp. Here "two component escape" refers to the com bined neut ron and direct escape components. The source redshift evolution is 
assumed to follow the GRB rate (star formati on rate from IHopkins Beacom (2006) including Kistler et al. (2009) correction for GRB 
evolution). The HiRes data are taken from lAbbasi et al.l ( 120081) . in this example it is assumed that all bursts are alike in the co- moving 
frame with r) = I, U = 3.3 ■ 10"^ s, F = lO^-S, L^^^^^ = 7 ■ 10^1 erg-s"!, e'^ ^^^^^ = 14.76 keV, a = 1, /3 = 2, and fc = 1 in equation 

where t^, is given in the co-moving frame. The resulting maximal proton energy with these parameters is 1.9 • lO^'^ GeV (co-moving frame). 
The normalization of the spectra is arbitrarily chosen, i.e., not based on a fit. 




The "transition (or ankle) model": using an injec- 
tion index at the source ap ~ 2.0. It can reproduce 
the CR spectrum above the ankle, i.e., the region 
above about 10^" GeV, but the energy range be- 



low the ankle requires the transition to a dilferent 
(galactic?) component. 

Normally, the escape mechanism from the source is not 
specified, only the cosmic ray injection spectrum (with 
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index ap and maximal proton energy) . We directly com- 
pute the cosmic ray direct escape and neutron escape 
spectra with our source model. Note that Up is the pro- 
ton injection spectral index in our model, not the cosmic 
ray injection index. Since multi-pion processes make the 
neutron spectrum somewhat harder (see, e.g., Figure S]), 
there is a slight difference compared to the usually as- 
sumed injection index. Our index ap refers directly to 
the acceleration processes in the source. 

Here we propose a third possibility: a two- 
component escape model for GRBs with both the 
neutron and direct escape components. The result is 
illustrated in Fig. [S] as a black thick solid curve. For 
ap = 2.5, the shape of the cosmic ray spectrum around 
dip and ankle is almost perfectly reproduced. It is obvi- 
ous that because of the large ap, the model works better 
than the transition model with ap = 2.0 (blue dotted). 
The red dashed curve represents a dip model with only 
the neutron escape {ap — 2.3), which obviously repro- 
duces the generic features of the observed spectrum, but 
not as pronounced as the two-component model. A com- 
parison of the two curves for ap = 2.5 (two component 
versus neutron escape only) illustrates the effect of the 
second escape component: it enhances the high energy 
part. The enhancement depends on the chosen parame- 
ter set and value of the acceleration efficiency. Note that 
Bethe-Heitler losses also contribute to the dip in the two- 
component model, but the effects of Bethe-Heitler losses 
and spectral break in the injection spectrum add up. 

We have also tested the impact on the cosmogenic neu- 
trino flux. A two-component model as shown in Figure [5] 
for ttp = 2.5 in combination with the GRB strong evolu- 
tion would lead to a very large productio n of cosinogeni c 
neutrinos close to the current bound (jlshiharal I2012D , 
which should be testable in the very near future. This 
is a generic feature common to all models with strong 
source evolution and a pronounced rise after the ankle in 
the cosmic ray spectrum; see, e.g., iKotera et al.l (|2010( ). 

8. SUMMARY AND CONCLUSIONS 

We have discussed the escape of UHECR protons from 
GRBs in the internal collision (prompt) phase in the 
framework of the GRB fireball model. We have iden- 
tifled three different regimes in terms of the dominating 
UHECR injection mechanism: 

1. Optically thin to neutron escape regime. 

Neutrons from photohadronic interactions, which 
are not magnetically confined, can escape from the 
source ( "neutron model" ) . 

2. Direct escape regime. Directly escaping protons 
from the outer edges of the shells dominate the 
UHECR injection, at least at the highest energies. 

3. Optically thick to neutron escape regime. 

Only neutrons from the outer edges of the shells 
can escape. 

In case (1), one (muon) neutrino per cosmic ray will be 
produced; in case (2), the UHECR escape will not be 
necessarily accompanied by neutrino production; and in 
case (3), the neutrino production will be enhanced com- 
pared to case (1), since the neutrinos can escape from 



everywhere within the shell. Normally, case (1) has been 
assumed in the literature. 

As one of our main results, we have identified the GRB 
parameter regions which are associated with these three 
different cases. We have demonstrated that direct es- 
cape is, for instance, important for GRBs with large F, 
such as we have shown for several typical ferTni-LAT 
observed bursts. The actual magnitude of the direct es- 
cape component somewhat depends on the acceleration 
efficiency, since protons can only escape directly at the 
highest energies, where the Larmor radius is of the order 
of the shell width. For efficient proton acceleration, we 
have demonstrated that the standard case (1), one neu- 
trino per cosmic ray, only applies to a very narrow region 
of the parameter space at the highest energies, since ei- 
ther direct escape dominates, or the optical thickness to 
neutron escape is large. For less efficient proton accelera- 
tion, a significant region where the standard assumption 
(1) applies has been found, which is around the often 
assumed standard parameter values. However, we have 
illustrated that the maximal proton energies are in that 
case not sufficient to describe the observed UHECR spec- 
trum for typical burst parameters. Therefore, it appears 
that the standard case (1) is in tension with the assump- 
tion that GRBs are the sources of the UHECRs, and spe- 
cific conclusions can be only drawn on a burst-by-burst 
basis. Note, however, that even in the direct-escape dom- 
inated cases, the neutron escape will contribute at lower 
energies as well, which leads to a spectral break in the 
cosmic ray injection spectrum. 

The region in which the optical thickness to neutron 
escape is large, case (3), has been found to be larger 
than previous calculations suggest, because high-energy 
processes have been included in the photohadronic in- 
teraction rate. In this case, the neutrinos can escape 
from everywhere within the shell, whereas the neutrons 
(and protons) are trapped over the photohadronic inter- 
action length scale and can only (directly) escape from 
the edges. The neutrino production can therefore be 
significantly enhanced. Furthermore, we have explicitly 
demonstrated that the formula used for the pion produc- 
tion e f ficien cy in the IceCube treatment in lAbbasi et al.l 
()2010l I2012D in fact underestimates the neutrino produc- 
tion in t he optically thick ca se, and that the original for- 
mula in iGuetta et al.l (|2004D applies instead (which was 
originally developed for the optically thin case) . This has 
consequences for individual GRBs, and is, in fact, already 
taken into account in the prediction by IHiimmer et al.l 
(|20T1) . 

From the model perspective, we have described the 
direct escape of protons and neutrons in a unified frame- 
work, and we have shown how the expansion of the fire- 
ball affects the escape. For instance, we have demon- 
strated that for a relativistic plasma with an adiabatic 
index 7 = 4/3 the ratio between shell width and direct 
escape mean free path is independent of the shell ra- 
dius, which means that the conclusions do not depend 
on the details of the time evolution. Note that there 
may be additional escape components, such as from dif- 
fusion, which however require a time-dependent calcula- 
tion, whereas the direct escape component is guaranteed 
in a rather model-independent way. Such additional com- 
ponents would however not affect our main arguments. 
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As far as the consequences for the UHECR observa- 
tions are concerned, we have demonstrated that a two- 
component escape model, which includes cosmic rays 
from direct proton and neutron escape, can reproduce 
the regions around the dip and ankle extremely well, be- 
cause the transition between the two components leads 
to a spectral break in the cosmic ray spectrum from a 
single source. Since an additional diffusive escape com- 
ponent may affect this, and the effect depends on the 
GRB parameter set, it remains to be seen if such an ex- 
tension of the dip model survives in a parameter space 
study. We expect that future neutrino observations will 
provide stringent constraints on two-component models 
using both the source (PeV) and cosmogenic (EeV) neu- 
trino fluxes. While the source neutrino flux is correlated 
with the neutron production, the cosmogenic neutrino 
flux is insensitive to the escape mechanism of the cosmic 
rays; a detailed discussion of this issue will be performed 
elsewhere. 

Finally, we note that the three different encountered 
regimes may even be present in one source, especially 
if collisions occur at very different radii. Therefore, the 
one neutrino per cosmic ray assumption is, in fact, not as 
general as one may believe. Note that some of our con- 
clusions can be transferred to other classes of sources. 



such as active galactic nuclei, and to heavier nuclei ac- 
celerated in the sources. In all those substantial 
fraction of particles may directly leak from the sources 
at the highest energies, and there can be regions where 
the source is optically thick to baryon escape. 
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APPENDIX 
A. ON THE EFFECTS OF DIFFUSION 

The description of diffusion depends on the actual magnetic field configuration, relevant scattering processes, and 
properties of the plasma. This means that a model-independent approach is not possible and a dedicated treatment of 
the transport equations, such as a set of Fokker-Planck equations, is required. However, we make some estimates for 
the additional effects of diffusion in this appendix, and demonstrate that they do not affect our qualitative conclusions. 

The diffusion length A^jfj, over which the particles can escape, is given by 



Dt' 



dyn ■ 



(Al) 



on the dynamical timescale ^^yn- Here D is the diffusion coefficient D = L'^/T', where L' is the displacement scale and 
T' time between the collisions relevant for diffusion. Since acceleration is assumed to be present over t'^y^, energy losses 
will take over thereafter, and conclusions on the additional escape after t'^y,^ can only come from a time-dependent 
calculation. 

Now, the diffusion coefficient depends on the actual field configuration. For example, if the magnetic field is aligned on 
scales larger than the Larmor radius, photohadronic and other scattering processes can lead to diffusion perpendicular 
to the magnetic field. If the scattering by photohadronic interactions dominates and energy losses can be neglected, 
D ~ R'l/t' and therefore 



•^diff ■ 



(A2) 



PI 



at the highest energies. This corresponds to the direct escape component, enhanced by the optical thickness y/T\l 
(cf., Section U), since the photohadronic scattering processes lead to the transport of additional protons closer to the 
shell boundaries. From Figure [21 however, one can easily see that this additional enhancement does not affect the 
qualitative conclusion that neutron escape dominates in the optically thick regime. 

Another regime which may be relevant for the diffusive escape in GRBs is the case of Bohm-like diffusion in magnetic 
fields. Assuming that our particles move with c, and re- writing the result in terms of i?^, we obtain D 
In order to study the impact of such a component, we assume that equation (|2ip is changed into 



i?^c foe B'-^] 



y^,^{E') = min [Ar', ^ R'^{E') ct'^y^,ct;^{E')] , 



(A3) 



where equation (jAip has been used to compute the fraction of diffusively escaping protons. This form implies that at 
most all of the protons can escape, and that this diffusive component will be suppressed if photohadronic interactions 
lead to significant energy losses or scattering. Note that other energy losses and time-dependent effects have been 
neglected here, which means that this can only serve as a rough estimate. In addition, it is implied that the amount of 
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Fig. 7. — Cosmic ray flux components as in Figure[T]and Figure|2] including a (model-dependent) estimate for a diffusive escape component 
of protons (orange curves); see main text for details. 



electrostatic turbulence in the plasma is only moderate. For a higher turbulence, actual MHD calculations are needed. 
This means that corrections and even a factor-of-a-few deviations are quite likely, and even the energy dependence 
may change if that assumption does not hold. 

We show in Figure[7]the impact of the additional diffusive escape component, given our assumptions, for our standard 
examples. For the optically thin sources (to neutron escape) in the upper row, the diffusive component may be as 
large as the direct escape at the highest energies, i.e., almost all protons can escape. This is expected, since direct 
escape and diffusion describe the same limit there and the Larmor radius is comparable to the dynamical timescale at 



the highest energies. However, since X'^-^g i 



the energy-dependence is different from direct escape, which 



means that at lower energies the diffusive escape may actually dominate. For the optically thick sources in the lower 
row, the diffusive escape can be enhanced compared to the direct escape because t'^^^ > R'^ at the highest energies 

(c/.. Figured! lower panels, where the acceleration rate, corresponding to cR'j^^, is larger than t'^\). The spectral 
breaks come from the points where photohadronic interactions take over. However, the neutron escape component 
still dominates above 10^ GeV. 

This example illustrates that an additional diffusive escape component does not affect our conclusions, since (a) it 
does not change the results significantly at the highest energies and (b) it can only enhance the ejection of protons 
from the sources, which means that the region where the one neutrino per cosmic ray connection holds will be further 
reduced. However, it also demonstrates that our direct escape component can only be regarded as a lower estimate 
for the actual cosmic ray leakage from the sources. As far as the impact on cosmic ray observations in Section [7] is 
concerned, the qualitative conclusions (additional component can better reproduce dip model) remain unchanged, but 
the actual energy dependence may depend on the interplay of different escape components. Since the estimates in this 
appendix are quite model-dependent, we do not mix them with the main text. 
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B. PION PRODUCTION EFFICIENCY IN THE OPTICALLY THICK REGIME 

There have been several shght l y diff erent definitions of in the literature; see, e.g., (jWaxman fc Bahcalll 119971 : 
IGuetta et al.l 1200^ lAbbasi et alj feoiOl ). In words, is defined as the fraction of proton energy going into pion 
production due to photohadronic interactions. For a single particle, this fraction has to be /,r < 1, whereas for a 
density of particles, this is no longer obvious and depends on the exact definition. One self-consistent possibility is to 
use the injection rates Q' for the calculation, leading to the estimate 



''PI 

Q'p ~ Q'p 

Here, we use the fact that the pions are produced from the protons by photohadronic interactions at a rate t'p^^ with a 
pion inheriting the average amount of energy (a;^) ~ 0.2 from the parent proton. In the optically thin case t'^]^ > t'^^^ , 
the particle injection counters the escape over the dynamical timescale in the steady state, and we get Q'p = • iV^ 
from equation ([7]). When we insert this for Q'^ into equation (jBip we obtain 

fi%n = KM<M<^ , (B2) 

which is exactly the result used by IGuetta et al.l (|2004D . 

However, this is no longer valid in the optically thick case t'^]^ < t'j^^ . As one of the subtleties, one may want to 
take into account multiple interactions. If the neutron decay rate can be neglected, as it is typically the case, protons 
and neutrons can be treated as one baryon system with similar interaction rates for all baryons. In that case, the 
particles are confined in this baryon system, and photohadronic energy losses dominate. One can easily solve the 
kinetic equation describing continuous energy losses for a power law spectrum, in order to obtain Q'^ = t^^^ • Np ■ (a;^). 
By inserting this for Q'p in equation (jBl|) . the result for /^r simplifies to 

/Shick^l . (B3) 

i.e., the baryons lose all their energy into pion production. Hence, when using the injection rates Q' for the definition 
of /tt, no values larger th an one can be obtained and conservation of energy is shown. The definition from Guetta 
et al. (jGuetta et al.l 1200^ at least correctly describes the optically thin case, where injection rates are used for the 



definition of /^r. 

On the other ha nd, a second possibi lity to define f-^ is implied by the IceCube method to relate the neutrino and 
photon fiuences in lAbbasi et al.l (j201Cl| ). used for all current state-of-the-art stacking analyses. Here, the fluences are 
related to each other by energy partition arguments, which correspond to energy p artition of the part icle densities in 
the source. The pion production efficiency is calculated in the conventional way (jGuetta et al.ll2003 ). modified such 
that /tt < 1 in the optically thick case. However, the definition of /^r implied in this approach is based on the particle 
(steady state) densities N': 

fi'^^§ ■ m 

For N^, we can estimate that « Q'^-t'^^^ as long as the dynamic escape rate is larger than the decay rate, 
^dyn > *dec' large cuough energies. For the sake of simplicity, we neglect synchrotron losses here, which of 

course leads t o an a dditional correctio n of the pion productio n efficiency, as discussed earlier in lHiimmer et al.l (j2012l ) 
as well as|L| (j2012b[ ) (included in cs in lHiimmer et al.1 (j2012f) ). This however does not change the argument here, but 
is included in our numerical calculations. In the optically thin case, we again have Np ~ Q'p- t'^y^, leading to 

„(2) ^ Q'tt ■ t'dyn _ % 

./ IT, thin ~ r)l il ~ ni ' ^ 
' ''dyn ^p 

which is equivalent to the definition in equation (jBll) and leads to the result in equation (jB2p . It therefore does not 
make a difference what approach is used. However, in the optically thick case, we again have N'p k. Q'p ■ t'p^ ■ (xt^)^^ 
(for the baryon system), and we obtain 

(2) ^ , '-^^ 'dyn I \ ^7 / \ _ ^Pl I \ /Bg^ 

Wp ■ l-p-y Hp tdyi-, tdyn 

where we have used in the last step that Q'^/Q'p ~ 1 in the optically thick case from above. Clearly, since t'^^^ ^ ^'dya 

is possible, /i-^'' can be larger than one. It is interesting that the form in equation (jB6j) is the same as in IGuetta et al.l 
(|2004,) . which means that this formula applies in the optically thin and thick cases, if the definition in equation (jB4p is 
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usedQ Since this applies to the IceCube treatment, the formula for /^r in lAbbasi et alj ()2010D can be easily adjusted to 
describe the optically thick regime correctly. The large neutrino flux in Figure [2] reflects the previously underestimated 
neutrino production in the opticall y thick case, whereas neutrons are trapped by photohadronic interactions. Note 
that numerical predictions, such as iHiimmer et all (|2012[ ). include this enhancement automatically, but still find an 
order of magnitude lower predictions than with the original method. The reason is that not all burst s are affected 
by th is. Especially two of the five most dominant bursts (of the original method) of the IC-40 sample (jAbbasi et al.l 
120111 ) are severely affected by this, which changes how these bursts contribute to the total fiux. These two bursts have 
a neutrino fiux which is enhanced by a factor of more than four in our approach. Another burst is mildly enhanced, 
i.e., by about 10%, while the other two originally dominating bursts are unaffected by this. As a result, only three of 
the original five bursts dominate in the new analysis. 

Now one may ask: where does this qualitative difference come from, and why does that not violate energy con- 
servation? If the energy partition arguments are applied to the actual densities in the source (and not the injection 
spectra), as it may be the more plausible approach at least for electrons and photons, the secondary production can 
become very large since it is proportional to the product of proton and photon densities; c/., equation The proton 
density in the source will be balanced by the injection Q'^ = i^"^ • ■ (a;^) in the steady state in the optically thick 
case. That is, for fixed (from energy partition), it is implied that more protons have to be injected if t^"^ increases. 
Thus, the definition in equation (jB4|) is with respect to an open box to which energy will be injected from the outside, 
whereas the definition in equation (|B1[) includes energy conservation automatically. Since in the optically thick case 
the maximal proton energy is typically limited by the photohadronic interaction rate, the proton acceleration will 
always be more efficient than the energy losses up to the maximal proton energy. This situation can be expected in 
the collision phase. If, however, the acceleration is switched off at some point, the proton density will decay quickly 
and energy conservation applies again. Note that by using the actual proton and photon densities as input for the 
photohadronic interactions, we automatically have the optically thick case included if we define energy equipartition 
with respect to the baryon system (i.e., do not distinguish between protons and neutrons in the source). In practice, 
there will be about 50% of protons and 50% of neutrons in a source with a large optical thickness, according to the 
secondary multiplicities of protons and neutrons. For neutrinos from optically thick sources, only observables sensitive 
to the tt'^/tt" ratio will be sensitive to the actual distribution between protons and neutrons, such as the Glashow 
resonance; see discussion in IHiimmer et al.1 (|2010a| ). 
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